Line Broadening and Decoherence of Electron Spins in Phosphorus-Doped Silicon Due 

to Environmental 29 Si Nuclear Spins 
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Phosphorus-doped silicon single crystals with 0.19 % < / < 99.2 %, where / is the concentration 
of 29 Si isotopes, are measured at 8 K using a pulsed electron spin resonance technique, thereby the 
effect of environmental 29 Si nuclear spins on the donor electron spin is systematically studied. The 
linewidth as a function of / shows a good agreement with theoretical analysis. We also report the 
phase memory time Tm of the donor electron spin dependent on both / and the crystal axis relative 
to the external magnetic field. 
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Electron spin resonance (ESR) in phosphorus-doped 
silicon (Si:P) at low temperatures is now a textbook ex- 
ample of how an electron spin in a solid interacts with 
nuclei. In this system, the donor 31 P nucleus, the elec- 
tron bound to it, and 29 Si nuclei randomly occupying 
/ = 4.7 % of the lattice sites of natural Si, all carry 
the smallest possible spin angular momentum h/2. The 
ESR spectrum exhibits a doublet separated by 4.2 mT 
due to the contact hyperfine (hf) interaction with 31 P. In 
addition, the electron is subject to the hf field from 29 Si 
nuclei, since the orbital part of the electron wavefunction 
spreads over thousands of lattice sites. It was soon real- 
ized that this interaction inhomogeneously broadens the 
ESR lines. A 29 Si nucleus at the Ith site shifts the preces- 
sion frequency of the electron spin by ±a;/2, according to 
the direction of the nuclear spin, ai is the isotropic hf con- 
stant at the Ith site. As long as the number of the electron 
spins participating in the ESR line is large and the total 
amount of the hf field each electron spin feels is stochas- 
tic, the resultant line shape will be approximated by the 
Gaussian distribution, for which the FWHM linewidth 
AB e and the root-mean-square linewidth A_B rms are re- 
lated by 



AB e = 2V2 In 2 AB, 



(1) 



AB rms is written, neglecting the polarization of 29 Si nu- 
clear spins, as 



AB rms — 



1/2 



(2) 



where a/ is rescaled in the unit of tesla, and the sum runs 
over all the sites but the origin. Electron-nuclear double 
resonance (ENDOR) has allowed the experimental deter- 
mination of ai at a number of sites, and a good agreement 
was found between the linewidth in natural Si and the 



one deduced from Eq. ((2J) Q. However, neither is obvi- 
ous, when / is reduced, if the square root dependence on 
/ still holds nor the line remains Gaussian. This is one 
of the subjects we will verify in this contribution. 

The inhomogeneous line broadening is a consequence 
of static shifts of the electron spin precession frequen- 
cies. In reality, 29 Si nuclei which caused the shift are 
mutually coupled through the magnetic dipolar interac- 
tions. Hence the nuclei are also able to cause a tempo- 
ral change of the electron spin precession frequency by 
flip- flopping their spins. This leads to decoherence of 
the electron spins due to spectral diffusion Q , which has 
gained revived attention in connection with spin-based 
quantum computing. It is suggested that decoherence 
due to nuclear-induced spectral diffusion could be an ul- 
timate bottleneck for quantum computing schemes based 
on GaAs quantum dots (QDs) 0,01 El- In this context, 
Si:P is a promising testbed for studying nuclear-induced 
spectral diffusion, not only because quantum computers 
based on Si:P have been proposed |7j, but also because 
the concentration of environmental nuclei is controllable, 
which is impossible in III-V compounds whose available 
isotopes all possess non-zero nuclear spins. 

In this contribution, we experimentally study inhomo- 
geneous line broadening and nuclear-induced decoherence 
of the P donor electron spins in isotopically controlled Si 
single crystals whose 29 Si concentrations / range from 
0.19 % to 99.2 %. Isotopically purified Si is now in wide 
use for the study of isotope effects in Si 8j , and prepa- 
rations of the samples with controlled / (such as 50 %) 
were accomplished in the following manner: First, two 
bars of isotopically purified 29 Si or 28 Si and natural Si 
were set together, and partially melted into one by opti- 
cal heating. The mass ratio of the two bars were chosen 
to achieve desired /. The combined bar was made single- 
crystal by the floating zone (Fz) method using another 
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TABLE I: Summary of isotopically controlled Si single crys- 
tals. The 29 Si concentration / [%], crystal growth method 
(Czochralski or floating zone), net donor concentration Nd 
[10 15 cm -3 ], spin-lattice relaxation time Ti [ms], FWHM of 
the ESR lines AB e [mT], and R = M 4 /3(M 2 ) 2 , obtained from 
experiments (R cxp ) and simulations (R s im), are listed. 



Sample 


r 


Cz/Fz 6 N d c 


Ti 


AB e 


-^cxp 




29 Si-0.2%« 


0.19 


Cz 


0.3 


11.0 


0.023 


1.9 


_! 


29 Si-0.3% 


0.31 


Fz 


0.6 


13.1 


0.082 


1.7 


_f 


29 Si-l% 


1.2 


Fz 


0.7 


10.0 


0.078 


1.6 


1.65 


29 Si-5% e 


4.7 


Cz 


0.8 


15.8 


0.26 


1.0 


1.15 


29 Si-10% 


10.3 


Fz 


1.6 


14.6 


0.42 


1.0 


1.06 


29 Si-50% 


47.9 


Fz 


0.6 


14.9 


0.89 


1.0 


1.00 


29 Si-100% dE 


99.2 


Cz 


0.8 


4.4 


1.22 


1.1 


0.98 



"Determined by secondary ion mass spectroscopy (SIMS). 

'Crystal growth by the Fz method was carried out with the same 
apparatus as used in Ref. 

c Determined from the temperature dependent Hall effect mea- 
sured by the van der Pauw method (down to about 20 K). Residual 
acceptor concentrations are of the orde r of 10 13 -10 14 cm~ 3 . 

Further information is found in Ref. [l3 29 Si-0.2% here corre- 
sponds to isotopically purified 30 Si in the reference. 

e The same samples as used in our previous work (Ref.HlTl. 

^The simulations exhibit significant counts at the origin (no nuclei 
are found in any of the 176 sites), and R is not calculated. 
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FIG. 1: (Color online) (a) Echo-detected ESR spectra of Si:P. 
Only the normalized high- field lines are shown. The hori- 
zontal (magnetic field) and vertical (echo intensity) axes are 
arbitrarily shifted, and the small background signals are sub- 
tracted for clarity, (b) FWHM linewidth AB e as a function 
of /. The dotted line is a theoretical prediction, as explained 
in the text. 



bar of natural Si as a seed crystal. All the procedures 
were carried out in high purity argon gas atmosphere. 
Care was taken to achieve uniform distribution of 29 Si 
atoms. Finally, the single crystals of about 5 mm in di- 
ameter and a few centimeters in length were grown along 
[001]. Information on the samples used in this work is 
summarized in Tabic [fl The net donor concentrations 
Nd were controlled to ~10 15 cm -3 to minimize unwanted 
dipolar interactions among the donor electron spins, but 
still ensuring the measurement sensitivity. 

Pulsed ESR experiments were carried out using a 
Bruker Elexsys E580 spectrometer equipped with an Ox- 
ford ER4118CF cryostat. We exclusively applied a Hahn 
echo sequence given by 7r/2-r-7r-T-Echo, where r is the 
interpulse delay and the duration of the 7r/2 pulse was set 
to 16 ns. The measurement temperature must be chosen 
properly 0]: It has to be below 10 K for the study of 
spectral diffusion, otherwise the phase memory time Tm 
is dominated by the spin-lattice relaxation time T\ due 
to an Orbach process 0, On the other hand, at 

temperatures below 4 K, extremely long T\ reduces the 
efficiency of our pulsed experiments in which each echo 
sequence needs to be repeated at a time interval much 
longer than T\ . We found 8 K to be optimal: low enough 
for Tm not to be limited by T\ but high enough to en- 
sure a reasonable measuring time. T\ listed in the fifth 
column of Tabic [IJ measured by an inversion recovery 
method, vary with samples, but show no marked correla- 
tions with / nor Nd- We presume that this difference in 
Ti has little effect on the difference in Tm, as long as T\ 
of one sample is much longer than its Tm- The repetition 



rate of the Hahn echo sequence was then set to 80 ms. 

The echo-detected ESR spectra were obtained by mea- 
suring the intensity of the Hahn echo as a function of 
the external magnetic field Bq with r fixed. The high- 
field lines of the doublets are shown in Fig. Q}(a). The 
lines are Gaussian in / > 4.7 %, while in / < 1.2 % the 
lines resembled Lorentzian (shown only for 29 Si-l%). The 
same features were confirmed from the low-field lines (not 
shown). Comparison between experimental lincwidths 
and those deduced from Eqs. and J2J is straightfor- 
ward. Hale and Mieher have determined a; for 176 lattice 
sites through ENDOR 01 > an d we will use the values re- 
ported in Table II of Ref. [pp. dividing them by g e fiBh. 
Figurenj-(b) plots both theoretical and experimental AB e 
as a function of /, which agree quite satisfactory. AB e at 
/ = 99.2 % is calculated to be 1.15 mT, which accounts 
for 94 % of the experimental value, 1.22 mT. The under- 
estimation of 6 % is acceptable, since we included only 
176 sites. The rest of thousands of sites will have very 
small isotropic hf constants, and increase the linewidth 
further, but not largely. This result, together with the 
fact that the linewidths are independent of the crystal 
orientations within experimental errors, justifies the ex- 
clusion of the anisotropic term of the hf interaction. 

The line shapes in /< 1.2 % are no longer Gaussian, 
and the use of Eq. is dubious. A general method to 
characterize a line shape is to calculate the quantity R 
defined as M 4 /3(Af 2 ) 2 , where M 2 and M 4 are the sec- 
ond and fourth moments of the line [J^. A Gaussian 
curve has R of unity. On the other hand, R for a pure 
Lorentzian curve cannot be defined, since both Mi and 
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FIG. 2: Simulations of the ESR lines of Si:P for / = 1.2 
%, 3.0 %, and 4.7 %. The total counts for each line is 10 5 , 
and the width of each bar is set to 0.005 mT. The horizontal 
(magnetic field) axes are arbitrarily shifted. 

M4 diverge. Hence R is a convenient measure to quantify 
the deviation from Gaussian. In the seventh column of 
Table R numerically calculated from the experimen- 
tal lines are given. We set moderate cut-offs to exclude 
the contribution from signals remote from the line cen- 
ters [l^. R deviates from unity for / < 1.2 %, and R — 

1.6 ~ 1.9 is consistent with the previous report by Feher 
in isotopically purified 28 Si of 99.88 % (R = 1.8 is given 
in Table V of Ref. 0). The reason for this change to hap- 
pen may be explained as follows: When / is reduced, 
it becomes rare to encounter 29 Si nuclei near the donor, 
where a; can be larger than the observed AB e 0] . Once 
such a 29 Si nucleus having a large a; is found, the elec- 
tron spin is kicked out from the center part of the line, 
and there is no chance of its returning to the center, since 
other 29 Si nuclei are likely to have small ai. Those elec- 
tron spins can only contribute to form wings of the line. 
This qualitative argument can be checked by a simple 
simulation, in which the known 176 sites are assumed to 
be occupied, independently with probability /, by 29 Si 
nuclei which are up or down with probability one-half. 
Then the hf shift of an electron spin is computed. Re- 
peating the procedure 10 5 times, a frequency histogram 
should portray the observed line. The eighth column of 
Table 0] lists R calculated from each simulated line, and 
Fig. 21 shows simulated lines for / = 1.2 %, 3.0 %, and 

4.7 %. It is seen that R gradually increases as reduc- 
ing /, and deviates from unity. The agreement between 
-R exp and i? s j m for / = 1.2 % is good, even though the 
histogram is strongly peaked at the origin due to the lim- 
ited number of the sites considered. It is interesting to 
continue this work to lower-/ samples in order to seek for 
the intrinsic linewidth, or the boundary at which other 
broadening mechanisms prevail fl7) ]. 

Now we examine the echo decay behavior by changing 
r. Bq was fixed at the centers of the high-field lines. We 
define Tm as the time at which the envelope of an echo 
decay curve damps to 1 /e of its initial echo intensity. Fig- 
ure [3] shows echo decay curves as a function of 2r, when 
Bq is along the [001] crystal axis. The strong oscillation 
observed in the high-/ samples is electron spin echo enve- 
lope modulation (ESEEM). Since this effect has already 
been analyzed by us [ll| and by Ferretti et al. [l8j , it will 
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FIG. 3: (Color) Hahn echo decay curves as a function of 2r 
when Bo is along the [001] crystal axis. Note that the scales 
of horizontal axes in upper and lower panels differ by an order 
of magnitude. Gaussian fits to the echo decay envelopes are 
superposed on the data of 29 Si-10%, 50%, and 100%. The 
oscillations due to ESEEM are also visible. 

not be dealt in the present contribution. 

For samples with / > 4.7 %, the echo decay envelopes 
are obscured by ESEEM. To deduce Tm, we approximate 
them as Gaussian of the form exp[— (2t/Tm) 2 ]- To deter- 
mine a more accurate form of the decay envelopes, which 
could contain several time constants of different origins, 
requires further research. We repeat the experiments as 
rotating the samples about the [110] axis to characterize 
the angular dependence of Tm- It takes maximum when 
Bq II [001] and minimum when Bq | [111] as shown in 
Fig. 0J- (a). This result bridges a large gap of / in our 
previous report of Tm carried out when only 29 Si-100% 
and 5% were available Tm also exhibits a power- 

law-like dependence on / as shown in Fig. @J (b). Empir- 
ical fittings suggest Tm oc / 7 with 7 = —0.86 for Bo || 
[001] and -0.89 for Bq || [111]. The 29 Si concentration 
dependence confirms that Tm is affected by the pres- 
ence of 29 Si nuclei, and the angular dependence demon- 
strates that spectral diffusion plays a major role, since 
the tendency reflects the strength of the nuclear dipolar 
couplings. When Bq | [111], one of the four tetrahe- 
dral bondings is parallel to Bo, and this pair of nuclei 
gives rise to the strong dipolar coupling, which makes 
Tm shortest. When Bq | [001], all the dipolar couplings 
between nearest neighbors vanish, so that Tm is longest. 

The angular dependence of T M in 29 Si-100% and 5% 
has been explicitly compared with theory developed by 
de Sousa and Das Sarma |4J in Fig. 1 of Ref. 19, show- 
ing reasonable agreement. There are other theories on 
nuclear-induced spectral diffusion, which are developed 
recently and applicable to Si:P as well as GaAs QDs [U|(J. 
We believe our experimental results serve as a testbed for 
the validations of such theories. 

The situation changes in / < 1.2 %. The forms of the 
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FIG. 4: (Color online) (a) Angular and (b) 29 Si concentration 
dependence of Tm- The circles (squares) in (b) represent Tm 
obtained when Bo is along [001] (approximately [111]). Only 
the data points in the filled symbols are used for the fittings. 



spectral diffusion at fixed temperature and donor con- 
centration. 

In conclusion, we prepared isotopically controlled Si 
single crystals, and investigated the widths of the ESR 
lines and Tm of the donor electron spin as a function 
of the 29 Si concentration / (0.19 % < f < 99.2 %) at 
8 K. We confirmed the square root dependence of the 
linewidth in / > 4.7 %, and also analyzed the deviation 
seen in / < 1.2 %. It is found that Tm is proportional 
to the power of / in / > 4.7 %, while further research is 
necessary in / < 1.2 %. We believe our results on Tm 
serve as a testbed for the validations of theories dealing 
with nuclear-induced spectral diffusion, and will provide 
insight into other systems such as electron spins confined 
in QD nanostructures. 
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echo decay curves are quite dissimilar to Gaussian. The 
fittings by a single time constant are difficult, and Tm 
as the 1/e decay points are plotted in Fig. 0}(b). T M 
of 29 Si-l% when Bo | [111] roughly coincides with the 
value extrapolated from the empirical power law. It must 
be noted that we encountered phase fluctuations of the 
quadrature-detected echo signals, occurring at 2t longer 
than about 500 /is, which significantly distorted the echo 
decay curves. This means, in effect, all the echoes are 
generated along one predetermined direction of the ro- 
tating frame, but each echo is detected along a random 
direction. A very similar instrumental problem has also 
been reported by Tyryshkin et al. when they measured 
purified 28 Si (/ - 0.08 %, N d = 10 15 cm" 3 ) 13]. Since 
our apparatus is of the same design as theirs, we suspect 
the two problems arise from the same origin J^jJ. Apart 
from this problem, it is not surprising that decoherence 
mechanisms other than spectral diffusion, such as the 
dipolar interactions among the electron spins and the T\ 
effect, both of which we have assumed to minimize, come 
into play when / is reduced 0, 0] . Indeed, if we ex- 
trapolate the power law to / < 1 %, Tm will soon be 
comparable with T\ at this temperature. Ultimately, the 
temperature and donor concentration dependence of Tm 
will have to be examined in detail to fully characterize 
the decoherence in this realm of /, which is beyond the 
scope of the present contribution. Our aim here is to in- 
vestigate the decoherence dominated by nuclear-induced 
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